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Abstract 

We study the cosmological evolution of isotropic matter on an infinitely thin coni- 
cal codimension-two brane-world. Our analysis is based on the boundary dynamics of 
a six-dimensional model in the presence of an induced gravity term on the brane and 
a Gauss-Bonnet term in the bulk. With the assumption that the bulk contains only a 
cosmological constant A^, we find that the isotropic evolution of the brane- universe 
imposes a tuned relation between the energy density and the brane equation of state. 
The evolution of the system has fixed points (attractors), which correspond to a fi- 
nal state of radiation for = and to de Sitter state for A^ > 0. Furthermore, 
considering anisotropic matter on the brane, the tuning of the parameters is lifted, 
and new regions of the parametric space are available for the cosmological evolution 
of the brane-uni verse. The analysis of the dynamics of the system shows that, the 
isotropic fixed points remain attractors of the system, and for values of A^ which 
give acceptable cosmological evolution of the equation of state, the line of isotropic 
tuning is a very weak attractor. The initial conditions, in this case, need to be fine 
tuned to have an evolution with acceptably small anisotropy. 
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1 Introduction 



Cosmology in theories with branes embedded in extra dimensions has been the subject 
of intense investigation during the last years. The most detailed analysis has been done 
for brane world models in five-dimensional space [1]. The effect of the extra dimension 
can modify the cosmological evolution, depending on the model, both at early and late 
times. Both modifications can be interesting phenomenologically. For example, the early 
time modification may give less fine tuned inflationary parameters [2] and the late time 
modifications may shed light to the recent acceleration of the universe [3] . 

Much less has been done, however, for the cosmology of theories in six or higher dimen- 
sions with branes of codimension greater than one. This is because, unlike the codimension 
one case, these branes exhibit bulk curvature singularities which are worse than 5-function 
singularities. They then need some regularization (introduction of brane thickness) which 
makes the study of cosmology on them rather complicated [4]. An alternative way to 
study cosmologies of branes of higher codimension would be to consider corrections to the 
gravitational action, such as an induced curvature term on the brane [5] and a Gauss- 
Bonnet term in the bulk [6], which allow the brane to have a mild singularity structure 
(see also [7]). These thin brane cosmologies would have the additional advantage that the 
internal structure of the brane docs not influence the macroscopic cosmological evolution. 

The case of codimension two is particularly interesting because these branes have a 
special property. Their vacuum energy (tension) does not curve their world-volume but 
only induces a deficit angle on the branes [8]. There have been several attempts to utilize 
this property in order to self tune the effective cosmological constant to zero and provide 
a solution to the cosmological constant problem [9] . The requirement of simple singularity 
structure in this case means that the branes are purely conical. If higher dimensional 
gravity is conventional {i.e., if there is only a higher dimensional Einstein term in the 
action), it has been found that cosmological evolution on the brane is not possible, because 
only pure tension is allowed on the brane [10]. The way out of this problem, keeping the 
structure of the singularity simple, is, as mentioned before, to assume more complicated 
higher dimensional gravity dynamics, as e.g., the introduction of a Gauss-Bonnet term in 
the bulk [11] and an induced gravity term on the brane [12]. 

The latter modifications, have as a consequence the existence of a brane Einstein equa- 
tion which is purely four dimensional [11, 12] with the mere addition of a deficit angle 
dependent cosmological term. In our previous work [12], we have shown that apart from 
this Einstein equation, the bulk Einstein equations evaluated on the brane provide a con- 
straint on the matter evolution on the brane. If only an induced gravity term is added, 
this constraint corresponds to a tuning between brane and bulk matter. If in addition a 
Gauss-Bonnet term is included in the bulk, this constraint can either be a matter tuning or 
a dynamical equation depending on the symmetries of the spacetime metric. This happens 
because the constraint involves the Riemann tensor of the induced four dimensional metric. 
If the induced metric is isotropic, the Riemann tensor can be expressed in terms of the 
Ricci tensor and the scalar curvature and, through the brane Einstein equation, will give a 
tuning between the matter of the bulk and the brane. If the four-dimensional metric is not 
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isotropic, the Riemann tensor is independent of the other two known curvature quantities 
and the constraint gives rise to a dynamical equation for the anisotropy. 

In the present paper, we will study in details the cosmological evolution of a conical 
brane with both an induced gravity and a Gauss-Bonnet term added in the higher dimen- 
sional gravity action. For simplicity we will assume that the only matter in the bulk is 
cosmological constant A^. We will solve the equations of motion evaluated on the brane 
and assume that the integration of them in the bulk does not give rise to pathologies {e.g., 
singularities). 

Firstly, we will study the isotropic cosmology, in which the brane matter has to obey a 
tuning relation. We will see that the evolution of the system for = tends to a fixed 
point with w — 1/3 and for A^ > to a fixed point with w — —1. For A^ < the system 
has a runaway behaviour to w — > oo. 

We will then relax the isotropy requirement for the metric (keeping, however, the matter 
distribution isotropic) in order to find whether the above matter tuning is an attractor or 
not. The matter on the brane need not now satisfy the previous tuning relation and the 
allowed regions of initial values of the energy density and pressure are significantly larger. 
The analysis of the dynamics of the system shows that line of isotropic tuning is an attractor 
for A^ > and thus the system isotropises towards it. However, for values of A^ which 
give a realistic cosmological evolution of the equation of state, the attractor property of 
the previous line is very weak and fine tuning of the initial conditions is necessary in order 
to have an evolution with acceptably small anisotropy. For Ab < the system shows, as 
in the isotropic case, a runaway behaviour. We will discuss in detail all the above cases 
and finally present our conclusions. 



2 Boundary Einstein Equations for a Conical Codi- 
mension-2 Brane-world Model 

We consider a six-dimensional theory with general bulk dynamics encoded in a Lagrangian 
jO-Buik and a 3-brane at some point r = of the two-dimensional internal space with general 
dynamics Cbrane in its world- volume. The gravitational dynamics is described by a Gauss- 
Bonnet term in the bulk and an induced four-dimensional curvature term localized at the 
position of the brane. Then the total action is written as 



/" 6 /" 6 

+ Id xCbuIU + d xCbrane ' 

where g^^pf is the six-dimensional metric, Mg is the six-dimensional Planck mass, M4 is 
the four-dimensional one, Tc — M^/Mq the cross-over scale between four-dimensional and 
six-dimensional gravity and o; > is the Gauss-Bonnet coupling constant. The singular 
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terms have been written in the particular coordinate system where the metric reads 



ds^ = g^,u{x, r)dx^dx'' + dr^ + L^{x, r)de'^ , (2) 

where gptv{x, 0) is the brane- world metric and denote four non-compact dimensions, ji — 
0, 3, whereas r, 9 denote the radial and angular coordinates of the two extra dimensions 
(the r direction may or may not be compact and the 9 coordinate ranges form to 27r). 
Capital M,N indices will take values in the six-dimensional space. Note, that we have 
assumed that there exists an azimuthal symmetry in the system, so that both the induced 
four- dimensional metric and the function L do not depend on 9. The normalization of the 
(5-function is the one discussed in [13]. 

The full equations of motion that are derived from the above action are 



rpiB)N rp{br)u;:t^ xA^ '^(^) 



(3) 



with G^^ ^ — R^^ ^ — \R^^^5^ the six-dimensional Einstein tensor, G^^^ — R^f^ ^ — 
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the four-dimensional Einstein tensor and 
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where the Gauss-Bonnet combination is 

Cqb = -R(6) - ^R^MN + ^MNKA ■ (5) 

In order that there are no curvature singularities more severe than conical, we will 
impose certain conditions on the value of the extrinsic curvature X^j, = gi^^ on the brane, 
where the prime denotes derivative with respect to r, and on the expansion coefficients of 
the function L 

L = pi{x)r + (32{xy + Psixy + . . . (6) 

These conditions read [11] 

K^,\r=0 = 0, (7) 

Pi — const. and P2 — . (8) 

Imposing these conditions and keeping only the finite part in L" /L, the Einstein equa- 
tions (3) can be evaluated at r = 0. The effective Einstein equations on the brane (obtained 
by equating the 5-function parts of the Einstein equations) are [11, 12] 

= [T^'' - A45a , (9) 
^^^pi 

with Mj,i = Ml[rl + 87r(l - (3i)a] and Aj = -27rM|(l - (3i). The various components of 
the bulk Einstein equations evaluated at r = are given in the following: 
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The {iiiy) component 




+4—r; + 2k'r; + k; 'r - 2k: 'r; - 2g,,K'^ 'r 



L" 



. (10) 



The (rr) component 



^-R - a^iR' - ARl^ + Rl^^,) = -^J^^> . 



(11) 



The {99) — (rr) component 




\ 1 (rp{B)e rp{B)r\ 



(12) 



The {/I, r) component 



(13) 



In the foUowing, we wiU study the above equations in a time dependent background, 
assuming that the bulk consists of a pure cosmological constant T^f = —Ab5^ and that 
the matter content of the brane is an isotropic fluid with J"^^^)*^ — S; diag(— pb, Pj,, Pb, Pb). 

3 The Constrained Isotropic Case 

Wc are interested in the cosmological evolution of a fiat isotropic brane-universe, therefore 
we will consider the following time dependent form of the metric (2) 



We can use the gauge freedom to fix A'"(f, 0) = 1, while we define A{t,0) = a{t). The curva- 
ture singularity avoidance condition (7) we imposed, dictates that A^'(t, 0) = A'{t,0) — 0, 
while the second derivatives of these metric functions are unconstrained. 
For this ansatz, the Einstein equations (9) and (11) give 



ds'^ = -N^t, r)dt^ + A^{t, r)dx^ + dr^ + L^ix, r)d9^ . 



(14) 




Mi ■ 



(15) 



(16) 



(17) 
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The equations (15) and (16) are the usual Priedmann and Raychaudhuri equations of 
a four-dimensional universe with a scale factor a, while the third equation (17) appears 
because of the presence of the bulk and acts as a constraint between the matter density 
and pressure on the brane. To see this, a simple manipulation of the above equations gives 



A, 



Mi 



1 2 A4 

2 ^3"m|; 



3Ph 



1 



2 A4 



PI 



+ 3"^l^£ — ' ^^^^ 



which shows a precise relation between A^, pi, and P5. To simplify the equations, we 
assume that the vacuum energy (tension) of the brane cancels the contribution A4 induced 
by the deficit angle, i.e., 



Pb = -A4 + Pm and = A4 + P„ 
with P^ = WcPm- Then the above equations read 

3^ = 



Ml, 



a a'^ 



-Wr 



while the constraint equation becomes 

Kb Pn 



Mi 



^(3wc - 1) + ^(3wc + 



PI 



(19) 

(20) 
(21) 

(22) 



Using the metric (14), the Einstein equations (10) and (12) give a set of bulk equations 
which involve the boundary values of the second r-derivatives of the metric at r = 



L" 
T 

~L 



N" 



1 + 12a- 



1 + 4q; 2- + - 



1 A" 


^ fa 




1 + 4a — + — 


.^'^ 


\ a a J _ 



Ag 
^2 Mi ' 



(23) 



N" 



1 + 4a ( h - 

a 

A I d d^ 

1 + 4a ( 

a 



a 

a 



r d^i 


A" 




1 + 12a^ 




l + 4a 2- + — 






\ a a'^ ) 



. 



(25) 



These equations can be solved for the second r-derivatives of the metric, as functions 
of the matter content on the brane, and in principle they can give us information about 
the structure of the bulk at r = 



T 



1 _|_ ArviiS 



[w, + 11 



Mi 



2a 



pm 



l + 2a^(we+l)(3we-l) 



(26) 
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L" 

T 



= 3 



A" 4««;c]^ 



- 1 



A Aa^ + 1 



M 



PI 



A, 



A" 



-]V^-3^(l + 2aK+l) 



(27) 
(28) 



A potential problem in the cosmology of the system would be, if the denominator of 



(26) is equal to zero, i.e., when Wc 



wf- with 



_ - (l + ± ^64a2^ + 32a^ + 13 



1 Orv-^ 



(29) 



When this happens, the six-dimensional curvature invariant will diverge close to the brane. 
Thus, after discussing the cosmological evolution of the brane world-volume, we should 
always check that it does not pass through a point which satisfies the above relation. 



4 Cosmological Evolution of an Isotropic Brane- 
Universe 

Prom the constraint relation (22), we can solve for Wc, the allowed equation of state of the 
matter on the brane. It should satisfy the following equation 

_9^S I Pm ( 1 _ 4 Pm \ 

= — ^ , (30) 

\ ^ ^i) 

with > so that the Hubble parameter is real. 

Before analyzing the system, let us note a first important difference between the system 
of the pure four dimensional dynamics and the one with the extra constraint added because 
of the presence of the bulk. In the four dimensional system, a constant w is allowed and 
its value is preserved during the evolution of the system. On the other hand, the evolution 
of the system with the extra constraint forbids any evolution with constant Wf.. Indeed, by 
differentiating (30) and using (20) and the conservation equation pm + 3(1 -|- Wc)Pra^ = 0, 
we can find a differential equation for Wc 

Wc + ^l + Wc)pm^ = 0. (31) 

Imposing a further condition of keeping Wc constant, would result to a constant p^, 
related to Wf. in a specific way, and by the conservation equation, to zero Hubble for a. 
Thus, an a priori fixing of Wc would result to an inconsistent system. 
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To study the cosmological evolution, we look at the system 

Pm = -3(1 + Wc)PmH , (33) 

where H = a/ais the Hubble parameter. We will analyze the above system of the isotropic 
case for As = and 7^ 0, because of the different features that arise in the two choices 
of this parameter. 



4.1 Evolution of the System for = 

From the above dynamical system, taking into account (30), we find that there is only one 
fixed point in the evolution, the one with 

(p^,i^^^/;e) = (0,0,l/3) . (34) 

Linear perturbation around this point reveals that it is an attractor. Prom (31) and the 
conservation equation (33) we find for the Hubble parameter 

« ^ f35) 

The above equation can easily be integrated and solved for Wc (the "— " sign in the solution 
of Wc is rejected because it gives imaginary Hubble parameter) 

= — + , . (36) 

3 3 VST^ ^ ' 

From this equation, we see that from any initial condition along the line of tuning (30), 
the expansion of the universe drives the equation of state to — >^ 1/3, i.e., radiation. 
We have also verified this by integrating the system numerically. During this cosmological 
evolution, it can never happen that Wc — wf (compare (29) with (30)) and thus the whole 
system is regular. 



4.2 Evolution of the System for 7^ 

Prom the dynamical system, as written in the previous subsection, taking into account (30) 
with As 7^ 0, we find that there is a fixed point in the cosmological evolution 

(p„,i/^w,)= (^p^,^,-l) , (37) 

with -^^^ = — f + ly^l + 1^^- Since we should have pm > 0, this fixed point exists only 
for As > and corresponds to a de Sitter vacuum. Linear perturbation around this point 
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Figure 1: The evolution of the equation of state Wc for the case of aA^/Mg = 10 with 
initial conditions Wc = .31 and = .02. 



reveals that it is an attractor. Thus, any matter density on the brane eventually evolves 
to a state of vacuum energy. During this evolution, the singular point Wc — wj can be 
encountered only if Ab/Mq > .576 (equating the expressions (29) and (30) the energy 
density be real and positive only for this range of A^). However, it is easily verified 

that for these "dangerous" values of A^, it is w~ < —1. Thus, even for those values of 
Ab for which the system is singular, the dangerous point w~ is reached only if initially 
Wc < < —1, i.e., for rather unphysical initial equations of state. All other evolutions 
will flow to the fixed point without ever passing though w~. 

A potentially interesting case, with a cosmological evolution resembling that of our 
universe, would be the one where < aAs/M^ <^ 1. Then, the line of isotropic tuning, 
as it is given by relation (30), has a maximum close to Wc ~ 1/3. The evolution of an 
initial energy density larger than the one corresponding to that maximum, will evolve 
towards Wc ~ 1/3, pass from Wc = and asymptote to Wc — > —1 (see Fig. 1 for an 
example). The asymptotic value of the effective cosmological constant at the fixed point 
will be oi-j^ ~ ^ ^ should be rather small to match with observations (when 
performing such a comparison, it is reasonable to assume that all the dimensionful scales of 
the theory are roughly of the same order, i.e., oT^I"^ ~ Mpi ~ Mg). The latter requirement 
of extremely small q;Ab/M| is the usual cosmological constant problem. Although the 
standard cosmological evolution is described by piecewise constant equations of state with 
w ~ 1/3, 0, —1, in the present theory the equation of state has always time dependence. 

For A^ < 0, there docs not exist any fixed point. The evolution of this system has 
a runaway behaviour and flows to — > oo while pm 0+. During this evolution, 
the singular point Wc — wf can be encountered only if Ab/Mq < —.326 (equating the 
expressions (29) and (30) the energy density Pm can be real and positive only for this 
range of A^). It is easily verifled that for these "dangerous" values of A^, it is > 1/3. 
Therefore even in the case in which A^ takes these values, the dangerous point is 
reached only if initially Wc < wf, 
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We finally note that, for Ab ^ 0, there does not exist a fixed point with (p^, if^) = (0, 0) 
because the equation of state Wc diverges at that point. 

5 The Unconstrained Anisotropic Case 

In the previous section we found some interesting cosmological evolutions of an isotropic 
brane-universe if > 0. We should keep in mind however, that in all cases studied, the 
energy density on the brane is tuned to the equation of state in a specific way, which seems 
at first sight artificial. Therefore, it would be worth studying some cosmological evolution, 
in a codimension-2 brane-world model, in which this tuning is not required. If the system 
then evolves towards the previously studied line of isotropic tuning, we will conclude that 
this tuning is an attractor and thus not artificial. 

To do so, we have to consider geometries which are not isotropic and in which the 
Riemann tensor cannot be expressed in terms of the Ricci tensor and the curvature scalar 
[12]. Then, the constraint equation (11) will not give rise to a brane-bulk matter tuning, 
but rather to a dynamical equation for the anisotropy of the space. For this purpose, let 
us consider the following anisotropic ansatz where the metric functions depend only on the 
time t and the radial coordinate r [i.e., keeping the azimuthal symmetry) 

3 

ds^ = -N'^{t, r)dt^ + J2 A'(^> {dx'f + dr^ + L\x, r)de^ . (38) 

1=1 

We can again use the gauge freedom to fix N{t,0) = 1, while we define Ai{t,0) = ai{t). 
The singularity conditions dictate as before that A^'(t,0) = A[{t,0) = 0, while the second 
derivatives of these metric functions are unconstrained. The most general anisotropic 
evolution scale factors with the above property can be written as 

Ai{t, r) = a{f)h{t)c{t) + a(t)r' + . . . , (39) 
^(^,r) = ^ + 6(^)r^ + -- - , (40) 

^3(t,r) = ^ + 6(t)r^ + ... (41) 

where a = (010203)^/^ represents the "mean" scale factor and 6, c represent two degrees of 
anisotropy. To simplify further the analysis of the system, we will choose c = const, (by 
a coordinate redefinition then we can always set c = 1). The dynamics of this particular 
choice can help us to understand the qualitative features of the general case. 

Note that, with the choice c = const., we have Ai(t, 0)yl2(t, 0) = ^K^, 0), but not also 
Ai{t,r)A2{t,r) = A^{t,r). The later relation would lead to extra constraints which will 
overdetermine the system. For this ansatz, the Einstein equations (9) and (11) give 

3^-^ - P'^^' (A2) 
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c? b b'^ 
b ¥ 



a 



b 

2- + - + T 



a 



a a 
2- + - 



+ 



+ 



^ab ^ -Pi, + A4 
ab ~ 

ab ^ -Pb + A4 
ab Ml 



a'- 



b^ 



'PI 
-Pb + A4 



a 



+ Aa 



a 



aa 
> 

a 



db^ a^b 

+ 2 2 

2 ab^ a^h^ 



- 2: 



b^a 



Ml 



(43) 
(44) 
(45) 
(46) 



The equations (42), (43), (44), (45) are the equations of four dimensional Einstein 

gravity with the previously postulated anisotropy, while the last equation (46) appears 
because of the presence of the extra dimensions and corresponds to the constraint equation 
of the isotropic case that we considered. It is easy to see that from the three equations 
(43), (44), (45), one is redundant. This happens because we have frozen the c degree of 
freedom. Keeping then two linear combinations of them, the brane Einstein equations take 
a simpler form 



2- + 4— 



a 
b^ 



db 



Pb + A4 
Pb-Pb + 2A4 



52 ~'~ ab 



= . 



(47) 
(48) 
(49) 



These are purely four-dimensional equations, while equation (46) coming from the extra 
dimensions is now dynamical, providing a Hubble equation for b 



62 



Pb + A4 

' 4M|.; 16q; 



'Mi 



+ 



+ 2 



Pb I 

Pb 



.1 _|_ Orv-^fe 



M 



PI 



(i 



+ 2q; 



PI 



(50) 

We will now assume, as in the previous sections, that the vacuum energy (tension) of 
the brane cancels the contribution A4 induced my the deficit angle, i.e., that 



Pb = -A4 + pm and Pft = A4 + , 
with Pm — wpm- After this simplification, the above equations read 

,2 52 



Pr, 



a" 



62 



a a 
2- + 4— 
a 

6 _ 62 ah 
6~62^ d) 



= {1-w) 



0, 



(51) 

(52) 
(53) 
(54) 
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while the equation coming from the extra dimensions becomes 



62 



AMI, 



+ 



32a 



fO^B Pm 

Ml + Ml, 



(3W - 1) + 2(2w + 1)q; 



= f{Pm,w) . (55) 



[The "— " sign in front of the square root has been rejected, since it always gives rise to 
imaginary Hubble for either a or b.] It is interesting to observe that the Hubble equation 
(52) for the "mean" scale factor a. after substitution of (55), has apart from the usual linear 
term in p (of the conventional four- dimensional cosmology), additional correction terms in 
p. This is similar to what happens also to five-dimensional brane- world models [1, 14] 
and is due to the presence of extra dimensions. This modification happens only in the 
anisotropic case. In the pure isotropic case the four-dimcnsional brane- universe feels the 
extra dimensions by only adjusting its energy density to its equation of state, but without 
any modification in the structure of the Fricdmann equation. 

The Einstein equations (10) and (12) give again a set of bulk equations which involve 
the boundary values of the second r-derivatives of the metric at r = 0. It is analogous to 
the system (23), (24), (25) and has solution identical with the isotropic case (26), (27), 
(28) with A" /Ai = A" /A. As noted there, we should make sure that any evolution of the 
system does not pass through the points with w — wf, where wf given by (29), since then 
the quantities A'-/Ai, N"/N, L" / L will diverse. In addition, because in the anisotropic 
case < is allowed, as it can be seen from equation (52), evolutions which pass from 
apm/Mli — —l/A will give singular N"/N and L"/L, and therefore should be avoided. 



6 Cosmological Evolution of an Anisotropic Brane- 
Universe 

Before analyzing the system, let us note again that in contrast to pure four dimensional 
anisotropic dynamics, where a constant w is allowed, in the present case, where there is 
an extra dynamical equation because of the presence of the bulk, w has to evolve. Indeed, 
by differentiating (55) and then using the four dimensional equations of motion and the 
conservation equation p^ + 3(1 -|- w)pm^ = we can find a differential equation for w 



df 

aw 



df 

2f-(l + w)p, ■' 



dpn 



- = . (56) 

a 



Imposing a further condition to keep w constant, would result to a constant pm related to 
w in a specific way and by the conservation equation, to zero Hubble for a. Thus, an a 
priori fixing of w would result to an inconsistent system. 

To have real Hubble parameters for a and b, pm and w have to lie in specific regions 
for which the following inequalities are satisfied 

/>0 , / + ^>0. (57, 
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Define the following boundaries of the allowed regions 



Wi 



W2 



-OilB _|_ 



(i - 1"^) 



Q Pm 



+ 



pm 



Pm 



M 



PI 



(58) 



(59) 



The quantity wi coincides with the line of isotropic tuning Wc given in (30). The inequalities 
(57) can then be re-expressed as conditions for pm and w: 

• For p„ > 0, we should have w > wi. 

• For < 0, we distinguish two cases 

- for -| < ^ < 0, we should have w < W2- 

- for < — I, we should have w > W2- 

To study the cosmological evolution, we look at the system 



1.. , ^ P 



H, = -SH^H,, (61) 

Pm = -3(1 + w)pmHa, , (62) 

where Ha = a/ a and Hf, = b/b are the Hubble parameters for a and b respectively. The 
third equation is the energy conservation equation in which only the Hubble parameter 
Ha for the "mean" scale factor a appears. Again we will analyze the anisotropic case for 
Ab = and 7^ and we will compare the cosmological evolution with the cosmological 
evolution of the tuned isotropic case. 



6.1 Evolution of the System with = 

There are three different regions in the (w,pm) plane where these inequalities are satisfied 
(sec Fig. 2c). This relative freedom to choose the matter on the branc should be compared 
to the tuning that happens in the isotropic case. Relaxing the isotropy condition, the 
system can have initial conditions in a vast region of the parametric space. The line of 
isotropic tuning is just the boundary of region I. Furthermore, in the anisotropic case we 
have two other allowed regions for p^ < where evolution is also possible. 

From the dynamical system (60), (61), (62) we see that there is only one fixed point 
and that it is the same with that of the isotropic evolution, i.e., 

{pm,HlHlw)^ {0,0,0,1/3) . (63) 
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Figure 2: The three regions of allowed pm and w for: (a) > |, (b) < < |, 
(c) ^j-jT = and (d) ^gf- < 0. The thick lines correspond to the tunings w = Wc oi the 
isotropic case. The thick dots represent the fixed point attractors of the evolution. The 
lines of long dashing denote the singular points w = wf and apm/Mpi = —1/4. 

Linear perturbation around this point reveals again that it is an attmctor. 

The presence of the previous attractor fixed point will drive the system towards a 
final isotropic state of radiation. The way in which this fixed point is approached from 
an arbitrary initial energy density, can tell us whether the line of isotropic tuning is an 
attractor or not. If the anisotropy monotonically decreases during the evolution, it means 
the the line of isotropic tuning is an attractor. 

In order to analyze the features of the anisotropic evolutions, we proceed numerically. 
We solve the system of the two second order equations (53), (54) and the two first order 
equations (56), (62) for the four functions a, 6, p^, w. The initial conditions for and 
w are such that they lie in the allowed regions of Fig. 2c and the initial conditions for 6, a 
are such that (52) and (55) are satisfied. We can check that the later two conditions are 
satisfied throughout the numerical evolution of the system, although imposed only once 
initially. 

To understand how the anisotropy involves we define the mean anisotropy by the fol- 
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Figure 3: The time evolution for = in the three regions for different initial conditions: 
Wo = 0, 1/3,2/3 for the region I, wq = —1 for the region II and wo = 3/2 for region III. 
The thick lines correspond to the boundaries (58), (59) of the three regions and the stars 
to the initial conditions. The thick dot represents the fixed point of the evolution. 



lowing quantity 



Our analysis shows that in region I, whatever the initial conditions are, the system 
slowly tracks the line of isotropic tuning and eventually evolves towards w ^ 1/3 (see case 
I in Fig. 3). The anisotropy decreases during the evolution and thus the tuning between 
Pm and w in the isotropic case is an attractor. However, it is a rather weak attractor. The 
reason for this is that, as we find numerically, the anisotropy falls like A ~ t"^/^, while 
at the same time the energy density is redshifting much faster pm ~ This is also 
evident from the example of the case I in Fig. 3, where the tendency to approach the line 
of isotropic tuning is very weak. Asymptotically, when the fixed point is approached, we 
have y4 — 0, as expected. These cosmological evolutions give rise to a regular geometry at 
r = for initial conditions w < (see Fig. 2c). 




(64) 



where Hi = hi/ai (with defined after (38)) and (if) 



\Y.\=lH^ = a/a . 
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The time evolutions for regions II (see case II in Fig. 3) and III (see case III in Fig. 3) 
have on the other hand different characteristics. This is because there is no corresponding 
isotropic fixed points towards they could evolve. 

In region 11, the system evolves towards the boundary line, is reflected back and evolves 
asymptotically to pm — >^ 0~ at w — > — oo. The anisotropy increases as the boundary line 
is approached and decreases A \pl at the asymptotics. Let us note also that although 
region 11 is connected to the fixed point w = 1/3, we have seen no evohition towards this 
point. These cosmological evolutions give rise to a regular geometry at r = for initial 
conditions on the left of the lines of long dashing in the region II of Fig. 2c. 

In region III, the system evolves towards the boundary line as w — > oo and apm/Mpi 
—3/4. The anisotropy decreases but tends to a constant value A > -\/2. These cosmological 
evolutions give always rise to a regular geometry at r = 0. 

From the above analysis we conclude that, the relaxation of the tuning relation between 
w and pm, has as a consequence the appearance of new branches of brane world evolution 
(regions 11, 111), while in region I the system tends quickly to the isotropic fixed point 
attractor. Furthermore, we observe that the anisotropy in all cases decreases much more 
quickly than in the four dimensional case with the same initial conditions but without the 
extra dimensional constraint. 



6.2 Evolution of the System with 7^ 

Studying the asymptotics for wi and W2 from (58), (59) we find that for 7^ 0, there are 
three intervals of with different shape of the allowed regions in the {pm,w) plane. The 
parameter space for these distinct cases is shown in Figs. 2a,b,d. 

From the dynamical system (60), (61), (62) we see that there are two fixed points in 
general. The first one is the same with that of the isotropic evolution, i.e., 

{pm, Hi Hi w) = (^p/, 0, -1) , (65) 

with = — I + 1^1 + and it exists only for A^ > 0. Linear perturbation around 

this point reveals again that it is an attractor. 
The second fixed point that we find is 

{pm,HlHlw) = (^-Ab^,0,^,iJ . (66) 

and it exists only for A^ > 0. Linear perturbation around this point reveals that it is a 

repeler. 

Whenever the previous attractor fixed point exists, the system will be driven towards 
a final isotropic de Sitter state. The way in which this fixed point is approached from 
an arbitrary initial energy density, can tell us whether the line of isotropic tuning is an 
attractor or not. To analyze the features of the anisotropic evolutions, we will again proceed 
numerically. 
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For ^f- > |, the three allowed regions of the (pm, w) plane are connected (see Fig. 2a). 
However, we can see numerically that any evolutions with initial conditions in one of 
the three regions, stays in that region. In other words, there can be no crossing of the 
= — I or = lines. In region 1 the system evolves quickly towards the line of 
isotropic tuning and tracks it until the attractor fixed point of w = — 1 is reached. We can 
see numerically that anisotropy falls to zero like A ~ exp 



0{l)^^t and thus the 

line of isotropic tuning is attracting very strongly the evolution of the system towards it. 
In region II, the system evolves towards 0~ and w — > +oo. Its anisotropy decreases 

and in the asymptotics it tends to A — > -\/2- In region III the system evolves to w — > +oo 
and apm/Mpi — > —3/4. The anisotropy decreases, but tends to a constant value A > 

For < ^ < |, two of the allowed regions of the {pmiW) plane are connected (see 
Fig. 2b). We can again verify numerically that any evolutions with initial conditions in one 
of these two regions, stays in that region, so that there is no crossing of the |^ = line. 
In region I the system tracks the fine of isotropic tuning and evolves towards the attractor 
fixed point ol w = —I while it isotropises (A — > 0). The anisotropy falls to zero again like 

A ~ exp —0{1),J^^ t and thus, the fine of isotropic tuning is an attractor. However, 
the strength of the attractor can be very weak in the cases when <cl. In region II, 

6 

the system evolves towards Pm — > 0~ and w — > +oo. Its anisotropy decreases and in the 
asymptotics it tends to A ^ In region HI the system again evolves to w +oo and 
apm/Mpi —3/4. The anisotropy decreases, but tends to a constant value A > 

For < 0, all allowed regions of the {pm,w) plane are disconnected (sec Fig. 2d). In 
region I the system evolves towards pm — 0"*" and w — > +oo. Its anisotropy increases and 
tends in the asymptotics to A — > y/2. In region II, the system evolves towards pm — > 0~ 
and ty — > — oo in the same way as in the = case. Its anisotropy decreases and in the 
asymptotics it tends to A ^ y^. In region HI the system again evolves to -u; ^ +oo and 
apm/Mpi —3/4. The anisotropy decreases, but tends to a constant value A > \/2. 

The evolutions in the three regions resemble the ones plotted in Fig. 3 for the = 
case. In all three regions the cosmological evolutions give rise to a regular geometry at 
r — when the initial conditions are such that the system never crosses the lines of long 
dashing in Figs. 2a,b,d. For example, this happens for the evolutions in the region I for 
Ab > and wj < —1 < w < w'^. 

From the above we conclude that, the relaxation of the tuning between w and pm, has 
as a consequence the appearance of new branches of brane world evolution (regions II, HI), 
while in region I the system tends to the attractor fixed points, whenever they exist. For 
the cases when A^ > the lines of isotropic tuning are attractors, with strength depending 
on the value of A^. 

Furthermore, we observe that the anisotropy in all cases, apart the one of region I for 
Ab < 0, decreases much more quickly than in the four dimensional case, with the same 
initial conditions, but without the extra dimensional constraint. In region I for A^ < 0, 
the anisotropy increases and is larger than the one of the purely four dimensional case. 

Let us now examine again the interesting possibility of < aAs/M^ -C 1 with the 
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Figure 4: The evolution of the equation of state w and the anisotropy A for the case of 
ctA^/Mg = 10"*^ with initial condition wq = 2/3. Until the fixed point is approached 
{w ~ —1), we have significant anisotropy A ~ 0{1) — O(10~^). 



transition between w — 1/3 to w = and finally to w = —1. As it can be inferred from 
the previous discussion, the system tracks the line of isotropic tuning and evolves towards 
the attractor fixed point with w = —1 (see Fig. 4). However, due to the small value of 
aA^/Mg, the line of isotropic tuning is a very weak attractor. Most of the evolution is 
rather anisotropic with A ~ 0{1) — O{10~^) until the fixed point is approached, in which 
region it drops to zero. 

This large anisotropy makes the cosmological evolution phenomenologically problem- 
atic. In order that the anisotropy is acceptably small, the initial conditions for the energy 
density and the equation of state should be fine tuned to lie very close to the line of isotropic 
tuning initially. 

In conclusion, by analyzing the anisotropic dynamics of the system wc have seen that 
the lines of isotropic tuning are attractors for A^ > 0, with A^-dependent attracting 
strength. The most phenomenologically accepted evolutions with < aA^/Mg <C 1 do 
not isotropise quickly enough and thus need a fine tuning in order to evolve with acceptable 
anisotropy. 



7 Conclusions 

In this work, wc studied the cosmological dynamics of a conical codimension-2 brane- 
worlds. Wc considered a theory with a Gauss-Bonnet term in the six-dimensional bulk and 
an induced gravity term on the three-brane. For simplicity, we considered that the bulk 
matter consists only of a cosmological constant A^ but the brane matter is general and 
isotropic. We then analyzed in detail the Einstein equations evaluated on the boundary. 

We studied the system first for an isotropic metric ansatz. As was noted in [12], in 
the pure induced gravity dynamics there is a tuning between the matter allowed on the 
brane and in the bulk. This tuning, when the matter in the bulk is only a cosmological 
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constant, gives a precise relation between the matter density on the brane and its equation 
of state. If a Gauss-Bonnet term is added in the bulk, the constraint equation giving the 
previous tuning is modified by the addition of a Riemann squared term. However, since for 
isotropic evolutions the Riemann tensor can be expressed in terms of the Ricci tensor and 
the scalar curvature, the conclusion about the presence of the tuning remains the same as 
in the induced gravity case. 

We therefore considered the combination of the two effects for the isotropic metric 
ansatz and made the further simplifying assumption that the brane contains a tension 
contribution which gives a four- dimensional Einstein equation without a cosmological con- 
stant. The dynamics of the system of the four-dimensional equations plus the constraint 
equation of the extra dimensions, depend on the value of A^. If > the fixed point 
has w = —1 and corresponds to a de Sitter vacuum. If A^ = 0, the system has a fixed 
point with equation of state w = 1/3. Both of the previous fixed points are attractors. If 
Ag < the system has no fixed point and the evolutions exhibit a runaway behaviour to 
w — > -|-oo. 

The tight restriction on initial conditions on matter on the brane that the constraint 
equation imposes in the isotropic case, motivated us to look at the case where the initial 
ansatz is anisotropic. In this way, we could determine whether the above tuning corre- 
sponds to an attractor. For this purpose, we considered a particular kind of anisotropy. 
Since the Riemann tensor in this case cannot be expressed in terms of the Ricci tensor and 
the scalar curvature, the previous studied constraint becomes now a dynamical equation 
for the anisotropy of the brane- world volume. The matter on the brane and its equation 
of state, need not lie on the line of isotropic tuning, but can have values in vast regions of 
parametric space which was previously forbidden. There are always three distinct regions 
of parametric space, one of which has as a boundary the line of isotropic tuning. 

The dynamics of the system of the four-dimensional equations plus the extra dimen- 
sional equation, depend again on the value of A^. If the system starts its evolution in 
the region of parametric space which has as a boundary the line of isotropic tuning, it 
tracks the latter line and isotropises towards the attractor fixed point with w = — 1 for 
Ab > 0, or the one with w = 1/3 for A^ = 0. In the two other regions the system has a 
runaway behaviour w ^ oo and does not isotropise. For A^ < 0, the system has always 
a runaway behaviour. The important result of this analysis is that the line of isotropic 
tuning, coming from the constraints found in [12], is an attractor. However, for values of 
Ab which give acceptable cosmological evolutions, a fine tuning is unavoidable because of 
the weak strength of the above-mentioned attractor. 

The above conclusions rely on the specific assumptions that we have made in order 
to simplify the dynamics of the system and might be altered if we change them. Let us 
remind the reader once again what has been assumed in the above analysis which may 
affect our final conclusions. Firstly, we consider only conical branes and in order to study 
cosmology on them we consider certain corrections to the gravity action (induced curvature 
and Gauss-Bonnet term). Secondly, we assume that the bulk matter has only the form of 
a cosmological constant. Thirdly, in studying the cosmological dynamics we consider that 
the tension contribution of the brane matter cancels the one induced by the deficit angle. 
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Finally, there is an implicit assumption in this work that the complicated bulk equations 
of motion, when integrated away from the brane, are well behaved {e.g., no singularities). 

Understanding the dynamics of codimension-2 branes is interesting because of the im- 
provement they may make to the cosmological constant problem. There have been several 
models of self tuning relying on the properties of codimension-2 branes, where the vacuum 
energy of the brane does not curve the brane world- volume without fine tuning. It would 
be interesting to see what the cosmology of these models would be by repeating the above 
analysis for a bulk content which is provided by these models. 
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